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ABSTRACT 


NASA’s current effort for the thermal environmental control of the Space Station 
Freedom is directed towards the design, analysis, and development of an Active Thermal 
Control System (ATCS). A two-phase, flow-through condenser/radiator concept has been 
baselined, as a component of the ATCS, for the dissipation of space station thermal load 
into space. 

The proposed condenser concept is similar to the single-phase deployable 
photovoltaic radiator. Here the heat rejection occurs through direct condensation of ATCS 
working fluid (ammonia) in the small diameter radiator tubes. Analysis of the condensation 
process and design of condenser tubes are based on the available two-phase flow models 
for the prediction of flow regimes, heat transfer and pressure drops. The prediction 
formulas employ the existing empirical relationships of friction factor at gas-liquid 
interface. 

A recent preliminary analysis of condensation in the radiator tubes at Johnson Space 
Center revealed that annular flow covers over 90% of the tube length. This analysis uses 
prediction formulas for pressure drop which are, in general, based on interfacial friction 
factor model for steam-water flow in large tubes. The contributions made by capillary 
waves and transverse radius of curvature to the instability of interface are ignored in the 
current models. 

The prediction of the interfacial friction factor and the major uncertainties associated 
with it play a vital role in modelling the annular flow regime. The condensation process in 
the proposed ATCS radiator includes two-phase flow in small tubes. The capillary waves 
present in the annular flow regime contribute to the liquid-vapor interfacial friction and 
hence have a major influence on the pressure drop and heat transfer. We have reviewed the 
present models of micro condensation in two-phase annular flow regime. The majority of 
formulations, describing the annular condensation process, are based on either on 
homogeneous or one dimensional separated flow models. 

This study is an attempt to investigate the stability of interfacial waves in two-phase 
annular flow. It contains the formulation of a stability problem in cylindrical 
coordinates. The contribution of fluid viscosity, surface tension and transverse radius of 
curvature to the interfacial surface is included in this formulation. A solution is obtained 
for Kelvin-Helmholtz instability problem which can be used to determine the critical and 
"most dangerous" wavelengths for the interfacial wave. The annular condensation is 
examined in this study and the general equations describing the condensation process is 
included. 
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INTRODUCTION 


Two-phase condensing phenomena occur in a variety of industrial applications 
including condenser systems associated with nuclear and conventional power plants, space 
power generations, thermal energy conversion systems, vapor compression refrigeration 
cycles, and chemical processing. The ability to understand and to model the principal 
physical mechanism associated with the condensing flow processes is of considerable 
importance. 

Design and optimization of operating conditions of condensers under space 
environment conditions require qualitative information including pressure drop, void 
fraction distribution, quality distribution, flow stability. The behavior of condensation 
flow phenomena and the transition between flow regimes must be analyzed and understood 
before condenser systems may be designed for space applications. The approach adopted 
by the majority of the investigators has been to first focus on the changes that could occur 
in the flow regimes and their boundaries and also at the pressure drop at reduced-g 
environment. Rezkallah [1, 2] suggests that in two phase gas-liquid flow under 
microgravity conditions there exist a significant shift in the transition boundaries when 
compared with flow regime maps of the normal gravity. Under those conditions there is a 
pronounced influence of surface tension. 

The dynamic characteristics of two-phase condensing flow have received very little 
attention compared to their boiling flow processes. Commonly, the flow analysis of 
condensation process is based on existing prediction formulas for two-phase flow models 
[3], These prediction relations, in general, have been developed for two-phase, two- 
component systems (e.g. air-water or steam-water), which have often been obtained in 
rather large diameter tubes under earth's normal gravity field. The most widely used 
technique for determining two-phase flow regimes in tubes is that of Martinelli [4], which 
assumes that the pressure drop of the separate liquid and gas phases are equal, and makes 
no assumption regarding the specific flow configuration. 

Research on condensation heat transfer in two phase flow under microgravity 
conditions is very limited, and experimental data for heat transfer coefficients under those 
conditions are sparse. Our understanding of many important physical phenomena in two 
phase flow remains incomplete. The microgravity, high-vacuum, low-temperature 
environment of space cannot be duplicated on earth , and it is too costly to place more than 
a few large-scale tests in space. Therefore, the heat transfer coefficient correlations for 
single or two-phase flow, and the pressure drop multiplier relations for two-phase flow 
under low gravity conditions must be determined analytically. 

Surface tension forces play a significant role in determining flow regimes in small 
diameter tubes even at 1-g conditions. The presence of an interface between two fluid 
phases can influence the motion of fluids and the pressure drop when the interface has a 
finite curvature (transverse curvature) that is different from that at equilibrium [5]. In most 
studies of two-phase flow the effect of transverse curvature on the interface has been 
ignored and the interface has been assumed to be flat However, the pressure drop in two- 
phase flow is principally due to the influence of the interfacial roughness elements. Small 
capillary-type waves on the interface (ripples and disturbance waves) which act as surface 
roughness give rise to an increased pressure drop. Since gravitational forces at 1-g are 
nearly dominant in such flows, surface tension has been assumed insignificant, and 
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therefore given little or no consideration. But, under reduced gravity conditions there is a 
pronounced influence of surface tension on the flow condition. These waves vary in length 
and in amplitude. Some amplitudes are several times greater than the mean film thickness. 
In two-phase flow in small diameter tubes, the interface is curved and deformed as a result 
of the forces acting on it. Therefore, the shape of the interface must be determined as part 
of the solution of the problem of two-phase flow. 

The condensation process inside tubes may include several flow regimes. The most 
common flow regimes are annular flow (with liquid film distributed essentially uniformly 
along the tube wall), annular-stratified flow (with a preponderance of liquid flow along 
the tube bottom due to the gravitational action), wavy annular flow (in which the shearing 
action of higher vapor flow initiates the growth of large amplitude waves, occasionally 
reaching the top of the tube), slug flow (in which the regular and periodic bridging of the 
tube diameter results in separated slugs of vapor moving the liquid). 

In a space environment there is no gravitational force acting on the liquid film. 
Consequently a stratified flow regime is not expected to exist. The surface wetting 
characteristics typical of refrigerant fluids causes uniform distribution of the condensate 
film along the periphery of the tube. That is, the flow is annular type. In annular 
condensing flow, the liquid film remains relatively thin, and has increasing velocity in the 
direction of flow, so that at high qualities one expects good uniform heat transfer along the 
entire tube periphery. The vapor velocity decreases as the vapor mass flow rate decreases 
more rapidly than the vapor flow area. Eventually the vapor reaches a blunt interface where 
all the vapor is condensed. At this point the liquid velocity is quite low since it now 
occupies the entire flow channel. 

NASA's current effort for the thermal environmental control of the Space Station 
Freedom is directed towards the design, analysis, and development of an Active Thermal 
Control System (ATCS) to remove all the electrical loads, plus any biological or chemical 
generated heat loads [6]. A two-phase, flow-through condenser/radiator concept has been 
baselined, as a component of the ATCS, for the dissipation of the thermal load into space. 
Here the heat rejection occurs through direct condensation of the ATCS working fluid 
(ammonia) in the condenser/radiator tubes (each 10.08 ft. long with a nominal diameter of 
0.09 inches). Design requirements for the condenser tubes are: 98% quality ammonia at the 
inlet and a minimum of 10°F subcooling at the outlet of the condenser tubes with the 
constraint that the pressure drop in the condenser panel should not exceed 1.4 Psid. 

The pre-integrated truss heat rejection system consists of two radiator arrays: one 
array rejects heat from the 58°F ATCS loop and the other rejects heat from the 35°F 
starboard loops. Each array consists of three Orbital Replacement Unit (ORU) radiators 
connected in parallel. Each ORU contains eight radiator panels and has two parallel 
condensation paths. 

A preliminary analysis of condensation in the ATCS radiator tubes at Johnson 
Space Center [7] revealed that annular flow covers over 90% of the tube length. The 
analysis were based on modified prediction formulas in [3] which contains a compilation 
of the existing empirical relationships for two-phase flow models for the prediction of flow 
regimes, heat transfer, gas-liquid interfacial friction factor pressure drops. These 
analysis revealed that the fluid flow is predominantly laminar in the liquid film and 
turbulent in the vapor core. 
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OBJECTIVES 


The basic objective of this effort is to study the microgravity condensation process 
in small diameter tubes to support NASA's current mission on analysis of the two-phase 
flow-through condenser concept for the Space Station Freedom. The specific goals of this 
work to establish the forces that contribute to the friction factor at the liquid-vapor 
interface, to predict pressure drop and heat transfer, and to examine the interfacial 
instabilities contributing to the transition of flow from annular to slug flow regime. These 
will require a detailed analysis of the capillary waves at the liquid-vapor. 

This study will be limited to examination of condensation process in the annular 
flow regime of the condenser tube and the examination of the capillary waves at the 
liquid-vapor interface and their contribution to the interfacial friction factor. In the current 
investigation, the modeling of the capillary waves at the liquid-vapor interface is based on 
potential flow theory (inviscid fluid flow). A more detailed and complete analysis requires 
the consideration of shear forces at the liquid-vapor interface. 

We propose to use a separated two-phase flow model to study the condensing 
process in small diameter tubes. This requires solving the continuity, momentum and 
energy equations for each phase and using appropriate boundary conditions at the interface 
to couple the two phases together. Thus, die shape of the interface and proper formulation 
of tangential and normal stresses at the interface will be the focus of this research. 
Oscillating pressure difference component across the interface (due to surface tension and 
transverse curvature) which contributes to flow instability will be included in the proposed 
analysis. 

Most studies conducted on two-phase flow condensation assume that the liquid- 
vapor interface is flat and ignore the influence of interfacial waves on pressure gradient It 
is important to examine the behavior of the interfacial waves in two-phase flow and their 
contribution to pressure drop. The stability of interfacial waves is also important in 
determination of flow regimes in two-phase flow situations. The majority of investigations 
on interfacial wave behavior are based on inviscid fluid models. These works usually do 
not consider the effect of surface tension on interface. Most analysis of the effect of 
surface tension on the interface consider two-dimensional waves in Cartesian coordinate 
and hence do not include the contribution of the transverse radius of curvature; this is an 
important parameter in modeling two-phase flow in cylindrical pipes, especially in small 
diameter tubes. 

The analytical treatment of surface waves originated with Helmholtz [8], Kelvin 
[9], and Rayleigh [10]. Other researchers have attempted to extend these studies. For 
example, Taylor focused his investigation on inertial instability of ideal homogeneous 
incompressible fluid. The instability of the horizontal interface between two ideal 
incompressible fluids dates back to G.I. Taylor [11]. The modem extension of viscid 
flows has been reviewed by a number of researchers who have suggested ways to include 
interfacial surface tension and the fluid viscosity. In recent years few studies [12,13] have 
considered the contribution of fluid viscosity on the “ most susceptible” or “most 
dangerous” wavelength and the corresponding frequency. 
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Lienhard and Dhir [12,13] studied the viscous hydrodynamics instability of peak 
pool boiling problem. The formulation of their instability problem was based on a model in 
Cartesian coordinate, which considered a vapor jet flowing in a pool of liquid. However, 
in their analysis they included the contribution of transverse curvature of vapor jet to 
instability as a form of an additional oscillating pressure difference component across the 
interface. They suggested that the pressure transverse pressure effect at the liquid-vapor 
interface may be presented as: 


Apff = 



( 1 ) 


where a is the surface tension, Rj is the radius of the jet, and (z, t) is the height of the 
the propagating interfacial wave. 


ANALYSIS 

In this work we will attempt to formulate the instability of interfacial wave in 
annular two-phase flow. In general, formulation of the problem will include the fluid 
viscosity and interfacial pressure gradient resulting from the surface tension and 
transverse radius of curvature. We will solve this problem for inviscid flow and make 
recommendations on the approach to solve the problem for viscous fluid flow. 


Formulatio n of Stability Problem: 

In formulating an interfacial wave model of annular flow we take an approach 
similar to that taken by Lienhard and Dhir [12,13]. However, we will use cylindrical 
coordinates to model our instability problem. The assumptions used in modeling the 
instability problem are: 

1 ) the liquid and vapor are incompressible 

2) both fluids are Newtonian 

3) the base flow is steady, axisymeteric and one dimensional 

4) The slope of interface is small everywhere 

5) the nonlinear effects are negligible 

6) the gravitational forces on the flow are negligible. 

A schematic model of interfacial wave in two-phase annular flow is shown in Fig. 
1. In this model vapor flows in the core of the cylinder and the liquid film flows adjacent 
to the wall. For this configuration we imagine a sinusoidal wave form at the interface and 
attempt to determine the propagation speed of the wave. The disturbance wave is assumed 
to oscillate as 


T| = £ e ik < z - ct ) 


( 2 ) 




eik(z-ct) 


where £ is the wave amplitude, k = 2 7t/X is the wave number describing the disturbance 
periodicity in the axial direction, X is the wavelength c = k to is the velocity of 
propagating wave and © is the frequency of wave. The equation of interface is 


r, =ri + T) = Tj + e e ik < 2 * ct > 


(3) 


where ij is the mean position of the interface. The criteria for the wave to be stable requires 
c to be real. If c is imaginary then the wave is unstable. The case c/i < 0 represents a 
decaying wave and the situation c/i > 0 denotes a growing wave. 


Modeling the annular two-phase flow as in viscid and axisymetric, the continuity 
and momentum equations reduce to: 


3ui , i d , . n 

_ + i_( rv i)= o 


(4) 
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dt dz dr 


dvj 

3t 




_L §Ei 
Pi dz 

_L §Pl 
P i dr 


+ Vi 


+ Vi 


fd 2 Ui 

d 2 ^ 

laz 2 

9r 2 


&. 

d 2 ^ 


9r2 



+ L<tl.YL 

r 9r t 2 


+ G r 


(5) 
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For the convenience and the space economy we have used the subscript "i" to 
represent the governing equations of both the liquid and vapor with a single set of 
equations (substitution of g and f for i gives relations for the vapor core liquid film region, 
respectively). Unless otherwise specified, we will use this subscript to represent the 
physical variables for both fluids with a single expression. 

Now let’s assume that steady annular flow Uj (r) is disturbed slightly such that 


Ui =Ui(r)+ u’i(r,z,t) 

(7) 

Vi = v'i(r,z,t) 

(8) 


where the perturbed components u* and v’ are infinitesimally small compared with the 
average velocity of the base flow (note that in condensation problem Uj = Ui (r, z) and one 
dimensional flow assumptions might not give accurate results). The pressure can be 
presented as 


where 


Pi=Pi(z)+p'i(r,z,t) 


(9) 


p' = E p'(r)e ik < z - ct ) 


( 10 ) 


We would like to determine whether the flow is unstable relative to the postulated 
disturbance. Substituting equations (7) through (10) into the governing equations (5) 
through (6), eliminating the steady state solution, and linearizing the results (ie. neglecting 
the terms of second order in u’j and v’j ) the following equations can be written for the 
perturbed flow 


9u'i 

Ik 


+ f| <”'■> =0 


( 11 ) 


3 “'i + U; 3U '' 


dt 


dz 


, 9Uj i dp'j 

+ v’i 3 -i- = --A- -£1. +Vi 
or ft dz 


, /ffu’j 9 2 u'j l 9 u'A 

l \9z 2 3r 2 r dr] 




1 dp'i , v fffv'i , d 2 v'i , 1 dv'i y±\ 

ft dr 1 3z 2 5r 2 r 3r r 2 ) 


( 12 ) 

(13) 


These equations must satisfy appropriate boundary conditions and interfacial 
requirements. For the vapor core the boundary conditions are 


v’g=0 


du'g 

dr 


= 0 


at 

at 


r = 0 


r = 0 


(14) 

(15) 


11-8 




Figure 2. Pressure oscillation due to the transverse radius of the curvature of the 
interface. 


The boundary conditions for the liquid film flow are 

v’ f = 0 at r = R (16) 

u’f = 0 at r = R (17) 

The following conditions are the result of force balance at the liquid-vapor interface. At 
the interface the normal stress exerted by the vapor is equal to those of the liquid film. 


_ UV „ U 1 1 U II | 

p 8 -2^-J- + a — + — = Pr -2 W 


r= n 


This force balance includes the surface tension contribution of pressure oscillation 
due to the transverse radius of the curvature of the interface as shown in Fig. 2. Note that 
the last term on the right hand side of equation (18) is the local time average of this 
pressure oscillation. A balance of tangential forces due to shear stresses of the liquid and 
vapor on the interface will yield 
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( 19 ) 


( 20 ) 


Now we introduce the perturbation stream function as 
y = eF(r)e ik < z - c ‘) 


such that 




and v' = - — 
r az 


Substituting equations (21) through (22) and (10) into equations (12) and (13) we 
get the following relations for the equation of motion 

r 3r r 9r pi i k r 2 Qfi j3 dr r dr J ^ 23 ) 

Eliminating p'i by combining equations (23) and (24) the following relationship will 
result. 

Vi ( d*Fi 2 , 3 ^Fi 3 9Fj ( 2 k 2 3Fi 2 k 2 ^^- + k 4 f] 

1 * c \3r 4r 9r 3 r 2 dr 2 i 3 ^r r dr dr 2 J 05 ) 


Equation (25) represents the "Qrr- Sommerfeld” equation in cylindrical 
coordinates. Equations (23) through (25) are similar to those developed in [14]. In 
fact the left hand side of equations (23) through (25) are identical to the 
corresponding equations in [14]. However, [14] omits terms in the left hand side 
of those equations corresponding to equations (23) and (25). 

It is more practical and useful for analysis of the instability problem to 
obtain a solution for equation (25) in dimensionless form . Let us define the 
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following dimensionless variables for the problem (for convenience we make 
dimensionless symbols identical to those used for the original physical variables). 


z 

R ; 

r R; 

£ = £ 
R ; 

„ 11 
n = R ; 

II 

r ‘ = R 

-=b° 

a 

c 

n 

II 

II 

P'i = P 1 2 
p-.Uo 

11 

0 


Where Uo is an average fluid velocity (e.g„ the superficial vapor velocity). Substituting 
these dimensionless variables in equations (23) through (25) results in the following 
relations: 


lm.-A ljM= , 1 (l PF; l ^Fj 1 3Fi dFi) 

r 1 1 ; 3r r f 3r Pl i Re, (r gfi r 2 r 3 3r r 3r j 

_J f^Ei + 3 dFj , 2d 2 dFj 2 +a*?\ 

i aRei \ ft* r ^3 r 2 ^2 ^ dr r dr ^2 j 


(26) 

(27) 


(28) 


where a = kR is the dimensionless wave number and Rej is the Reynold number defined 
as 


Re^ 

v, (29) 

The non-dimensional form of the disturbance wave, the equation of interface and 
the perturbation stream function are expressed by the following relations: 


T] = £ gi o (z - c t) 

(30) 

r, =iT + Tj =^+e e ia ( z * ct l 

(31) 

y = e F(r) e ia ( z * c 9 

(32) 


To analyze the stability of the interfacial waves we must solve equation (27), using 
relevant boundary and interfacial conditions. Note that for inviscid flow or high Reynold 
numbers the left hand side of the equation vanishes and it reduces to 
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( 33 ) 


(Ui-c) 




+ 



.^Us 
1 dr 2 


= 0 


The contribution of the radial variation of the base velocity to the instability 
problem is represented by the last bracket in equation (33). For uniform base flow this 
term also vanishes and equation reduces to Helmholtz instability problem in cylindrical 
coordinates: 




o 


(34) 


The relevant boundary and interfacial condition for the "Orr- Sommerfeld" problem, 
equation (27), can be obtained from equations (14) through (20). They are 


_SFg_ _3 ^F l 
8 3r 


at r = 0 


(35) 
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at r = 1 


(36) 


2 9V* i „ &r\ , OTiR 

L 2 av ’g|[Pf| 


Re s to pg RUq to 2 2 Pg Ujj rf 

l Ps Reg dr J\pgJ 

at r= n (37) 


at r= n 

(38) 

v’i = + 

dt dz 

at r= ri 

(39) 


There is no complete solution available for equation (27 ) at this time. Most existing 
solutions are for the special cases of the "Orr- Sommerfeld" equation. Other solutions treat 
the "Orr- Sommerfeld" equation in Cartesian coordinates. The existing treatment of the 
instability problems are aimed at flow conditions at low Reynold number. For example, 
the hydrodynamic instability studies conducted by Lienhard and Dhir [12,13] are based on 
the formulation of problem in Cartesian coordinates. Their perturbation method solves the 

instability problem for a case where Re « 1 and the wave number a is close to one. 

Solution of Kelvin-Helmholtz Instability Equation: 

Our next attempt is to determine the most “dangerous wavelength” by solving 
equation (34). For the vapor core we use the boundary condition, equation (35), and the 
kinematic condition at the interface, equation (39). For the liquid him we employ the 
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boundary condition, equation (36) and the kinematic condition at the interface, equation 
(37). 


The general solution for equation (34) is 

Fj = Ci rli(ar) + C 2 r Ki(ar) (40) 


where Ii and Ki are modified first order Bessel functions. Applying the boundary and 
kinematic conditions we obtain a solution for the vapor flow, expressed as 


Fg = - (U g - c) 


rli(ar) 

.It(axD. 


(41) 


The solution for the liquid film is 

K 1 (ar)I 1 (a)-I 1 (ar)K 1 (a) 


Ff=-(U f -c) 


. (ii (a q))(Ki(a)) - (K^a ^(Ijfa) ) 
From equation (26) we get the following expression forp’i 


(42) 


(43) 


Substituting this relation into equation (37), the continuity condition of normal stress at the 
interface can be presented as 


-(Ug-c 


.3F P 


dr 


Pg RUq 


GT) R 

az 2 + 2 p g Uq if 


= - (Uf- c ) 



(44) 


Substituting equations (30), (41), and (42) and solving for c, equation (44) will result in a 
complex number, where 


c = c R + ci (45) 

where cr and q are the real and the imaginary part of the wave velocity . The stability 
condition requires that 


q = 0 (46) 

The critical wave number, etc, and wavelength, Xc = 2 7t/ etc, can be determined 
from this condition. The most dangerous wavelength can be evaluated using the criteria 


d(-iac) _ n 
3a 


(47) 
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CONCLUSIONS AND RECOMMENDATIONS 


The condensation process in annular flow has been discussed. The major forces 
influencing the pressure drop have been identified. It has been noted that under micro 
gravity conditions, the forces that influence the flow include inertial, viscous, and surface 
tension. 

The stability of capillary waves in small tubes has been formulated. The general 
formulation of stability considers both viscous and surface tension effects. In addition, 
pressure oscillation due to the transverse radius of curvature of the interface has been 
included in this study. Kelvin-Helmholtz's Instability Problem has been solved. 
Equations (44) through (47) provide the most dangerous wavelength. 

A more complete solution of equation (28) is needed to determine the viscous effect on the 
stability of the interface. Most models for evaluation of friction factor for condensation 
process in the annular flow regime are based on one-dimensional, two-phase flow in large 
tubes under normal gravity conditions. These models ignore the effect of capillary waves 
on pressure drop. These models require a friction factor multiplier based on an 
experimental value. All our attempts of modeling the condensation process resulted in a 
need for an experimental value of friction factor multiplier to determine the pressure drop 
due to interaction at the interface. 

There is need for further study of problems considered in this work. This includes 
more in depth solutions to the instability problem and attempts to solve the pressure drop 
problem in two-dimensional flow. A more detailed and complete analysis is required to 
determine the influence of shear forces on the instability of the interfacial waves. 

A complete analysis of the condensation problem requires solving the continuity, 
momentum and energy equations for each phase and using appropriate boundary conditions 
at the interface to couple die two phases. 
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